High-throughput characterization of mutations in genes that drive clonal evolution using multiplex adaptome capture sequencing.

Understanding how cells are likely to evolve can guide medical interventions and bioengineering efforts that must contend with unwanted mutations. The adaptome of a cell-the neighborhood of genetic changes that are most likely to drive adaptation in a given environment-can be mapped by tracking rare beneficial variants during the early stages of clonal evolution. We used multiplex adaptome capture sequencing (mAdCap-seq), a procedure that combines unique molecular identifiers and hybridization-based enrichment, to characterize mutations in eight Escherichia coli genes known to be under selection in a laboratory environment. We tracked 301 mutations at frequencies as low as 0.01% and inferred the fitness effects of 240 of these mutations. There were distinct molecular signatures of selection on protein structure and function for the three genes with the most beneficial mutations. Our results demonstrate how mAdCap-seq can be used to deeply profile a targeted portion of a cell's adaptome.

Global Transcriptional Response of Methylorubrum extorquens to Formaldehyde Stress Expands the Role of EfgA and Is Distinct from Antibiotic Translational Inhibition.

The potency and indiscriminate nature of formaldehyde reactivity upon biological molecules make it a universal stressor. However, some organisms such as Methylorubrum extorquens possess means to rapidly and effectively mitigate formaldehyde-induced damage. EfgA is a recently identified formaldehyde sensor predicted to halt translation in response to elevated formaldehyde as a means to protect cells. Herein, we investigate growth and changes in gene expression to understand how M. extorquens responds to formaldehyde with and without the EfgA-formaldehyde-mediated translational response, and how this mechanism compares to antibiotic-mediated translation inhibition. These distinct mechanisms of translation inhibition have notable differences: they each involve different specific players and in addition, formaldehyde also acts as a general, multi-target stressor and a potential carbon source. We present findings demonstrating that in addition to its characterized impact on translation, functional EfgA allows for a rapid and robust transcriptional response to formaldehyde and that removal of EfgA leads to heightened proteotoxic and genotoxic stress in the presence of increased formaldehyde levels. We also found that many downstream consequences of translation inhibition were shared by EfgA-formaldehyde- and kanamycin-mediated translation inhibition. Our work uncovered additional layers of regulatory control enacted by functional EfgA upon experiencing formaldehyde stress, and further demonstrated the importance this protein plays at both transcriptional and translational levels in this model methylotroph.

Evolution of satellite plasmids can prolong the maintenance of newly acquired accessory genes in bacteria.

Transmissible plasmids spread genes encoding antibiotic resistance and other traits to new bacterial species. Here we report that laboratory populations of Escherichia coli with a newly acquired IncQ plasmid often evolve 'satellite plasmids' with deletions of accessory genes and genes required for plasmid replication. Satellite plasmids are molecular parasites: their presence reduces the copy number of the full-length plasmid on which they rely for their continued replication. Cells with satellite plasmids gain an immediate fitness advantage from reducing burdensome expression of accessory genes. Yet, they maintain copies of these genes and the complete plasmid, which potentially enables them to benefit from and transmit the traits they encode in the future. Evolution of satellite plasmids is transient. Cells that entirely lose accessory gene function or plasmid mobility dominate in the long run. Satellite plasmids also evolve in Snodgrassella alvi colonizing the honey bee gut, suggesting that this mechanism may broadly contribute to the importance of IncQ plasmids as agents of bacterial gene transfer in nature.

Directed evolution of Escherichia coli with lower-than-natural plasmid mutation rates.

Unwanted evolution of designed DNA sequences limits metabolic and genome engineering efforts. Engineered functions that are burdensome to host cells and slow their replication are rapidly inactivated by mutations, and unplanned mutations with unpredictable effects often accumulate alongside designed changes in large-scale genome editing projects. We developed a directed evolution strategy, Periodic Reselection for Evolutionarily Reliable Variants (PResERV), to discover mutations that prolong the function of a burdensome DNA sequence in an engineered organism. Here, we used PResERV to isolate Escherichia coli cells that replicate ColE1-type plasmids with higher fidelity. We found mutations in DNA polymerase I and in RNase E that reduce plasmid mutation rates by 6- to 30-fold. The PResERV method implicitly selects to maintain the growth rate of host cells, and high plasmid copy numbers and gene expression levels are maintained in some of the evolved E. coli strains, indicating that it is possible to improve the genetic stability of cellular chassis without encountering trade-offs in other desirable performance characteristics. Utilizing these new antimutator E. coli and applying PResERV to other organisms in the future promises to prevent evolutionary failures and unpredictability to provide a more stable genetic foundation for synthetic biology.

Specificity of genome evolution in experimental populations of Escherichia coli evolved at different temperatures.

Isolated populations derived from a common ancestor are expected to diverge genetically and phenotypically as they adapt to different local environments. To examine this process, 30 populations of Escherichia coli were evolved for 2,000 generations, with six in each of five different thermal regimes: constant 20 °C, 32 °C, 37 °C, 42 °C, and daily alternations between 32 °C and 42 °C. Here, we sequenced the genomes of one endpoint clone from each population to test whether the history of adaptation in different thermal regimes was evident at the genomic level. The evolved strains had accumulated ∼5.3 mutations, on average, and exhibited distinct signatures of adaptation to the different environments. On average, two strains that evolved under the same regime exhibited ∼17% overlap in which genes were mutated, whereas pairs that evolved under different conditions shared only ∼4%. For example, all six strains evolved at 32 °C had mutations in nadR, whereas none of the other 24 strains did. However, a population evolved at 37 °C for an additional 18,000 generations eventually accumulated mutations in the signature genes strongly associated with adaptation to the other temperature regimes. Two mutations that arose in one temperature treatment tended to be beneficial when tested in the others, although less so than in the regime in which they evolved. These findings demonstrate that genomic signatures of adaptation can be highly specific, even with respect to subtle environmental differences, but that this imprint may become obscured over longer timescales as populations continue to change and adapt to the shared features of their environments.

Tempo and mode of genome evolution in a 50,000-generation experiment.

Adaptation by natural selection depends on the rates, effects and interactions of many mutations, making it difficult to determine what proportion of mutations in an evolving lineage are beneficial. Here we analysed 264 complete genomes from 12 Escherichia coli populations to characterize their dynamics over 50,000 generations. The populations that retained the ancestral mutation rate support a model in which most fixed mutations are beneficial, the fraction of beneficial mutations declines as fitness rises, and neutral mutations accumulate at a constant rate. We also compared these populations to mutation-accumulation lines evolved under a bottlenecking regime that minimizes selection. Nonsynonymous mutations, intergenic mutations, insertions and deletions are overrepresented in the long-term populations, further supporting the inference that most mutations that reached high frequency were favoured by selection. These results illuminate the shifting balance of forces that govern genome evolution in populations adapting to a new environment.

Identifying structural variation in haploid microbial genomes from short-read resequencing data using breseq.

BACKGROUND: Mutations that alter chromosomal structure play critical roles in evolution and disease, including in the origin of new lifestyles and pathogenic traits in microbes. Large-scale rearrangements in genomes are often mediated by recombination events involving new or existing copies of mobile genetic elements, recently duplicated genes, or other repetitive sequences. Most current software programs for predicting structural variation from short-read DNA resequencing data are intended primarily for use on human genomes. They typically disregard information in reads mapping to repeat sequences, and significant post-processing and manual examination of their output is often required to rule out false-positive predictions and precisely describe mutational events. RESULTS: We have implemented an algorithm for identifying structural variation from DNA resequencing data as part of the breseq computational pipeline for predicting mutations in haploid microbial genomes. Our method evaluates the support for new sequence junctions present in a clonal sample from split-read alignments to a reference genome, including matches to repeat sequences. Then, it uses a statistical model of read coverage evenness to accept or reject these predictions. Finally, breseq combines predictions of new junctions and deleted chromosomal regions to output biologically relevant descriptions of mutations and their effects on genes. We demonstrate the performance of breseq on simulated Escherichia coli genomes with deletions generating unique breakpoint sequences, new insertions of mobile genetic elements, and deletions mediated by mobile elements. Then, we reanalyze data from an E. coli K-12 mutation accumulation evolution experiment in which structural variation was not previously identified. Transposon insertions and large-scale chromosomal changes detected by breseq account for ~25% of spontaneous mutations in this strain. In all cases, we find that breseq is able to reliably predict structural variation with modest read-depth coverage of the reference genome (>40-fold). CONCLUSIONS: Using breseq to predict structural variation should be useful for studies of microbial epidemiology, experimental evolution, synthetic biology, and genetics when a reference genome for a closely related strain is available. In these cases, breseq can discover mutations that may be responsible for important or unintended changes in genomes that might otherwise go undetected.

Large chromosomal rearrangements during a long-term evolution experiment with Escherichia coli.

UNLABELLED: Large-scale rearrangements may be important in evolution because they can alter chromosome organization and gene expression in ways not possible through point mutations. In a long-term evolution experiment, twelve Escherichia coli populations have been propagated in a glucose-limited environment for over 25 years. We used whole-genome mapping (optical mapping) combined with genome sequencing and PCR analysis to identify the large-scale chromosomal rearrangements in clones from each population after 40,000 generations. A total of 110 rearrangement events were detected, including 82 deletions, 19 inversions, and 9 duplications, with lineages having between 5 and 20 events. In three populations, successive rearrangements impacted particular regions. In five populations, rearrangements affected over a third of the chromosome. Most rearrangements involved recombination between insertion sequence (IS) elements, illustrating their importance in mediating genome plasticity. Two lines of evidence suggest that at least some of these rearrangements conferred higher fitness. First, parallel changes were observed across the independent populations, with ~65% of the rearrangements affecting the same loci in at least two populations. For example, the ribose-utilization operon and the manB-cpsG region were deleted in 12 and 10 populations, respectively, suggesting positive selection, and this inference was previously confirmed for the former case. Second, optical maps from clones sampled over time from one population showed that most rearrangements occurred early in the experiment, when fitness was increasing most rapidly. However, some rearrangements likely occur at high frequency and may have simply hitchhiked to fixation. In any case, large-scale rearrangements clearly influenced genomic evolution in these populations. IMPORTANCE: Bacterial chromosomes are dynamic structures shaped by long histories of evolution. Among genomic changes, large-scale DNA rearrangements can have important effects on the presence, order, and expression of genes. Whole-genome sequencing that relies on short DNA reads cannot identify all large-scale rearrangements. Therefore, deciphering changes in the overall organization of genomes requires alternative methods, such as optical mapping. We analyzed the longest-running microbial evolution experiment (more than 25 years of evolution in the laboratory) by optical mapping, genome sequencing, and PCR analyses. We found multiple large genome rearrangements in all 12 independently evolving populations. In most cases, it is unclear whether these changes were beneficial themselves or, alternatively, hitchhiked to fixation with other beneficial mutations. In any case, many genome rearrangements accumulated over decades of evolution, providing these populations with genetic plasticity reminiscent of that observed in some pathogenic bacteria.

Identification of mutations in laboratory-evolved microbes from next-generation sequencing data using breseq.

Next-generation DNA sequencing (NGS) can be used to reconstruct eco-evolutionary population dynamics and to identify the genetic basis of adaptation in laboratory evolution experiments. Here, we describe how to run the open-source breseq computational pipeline to identify and annotate genetic differences found in whole-genome and whole-population NGS data from haploid microbes where a high-quality reference genome is available. These methods can also be used to analyze mutants isolated in genetic screens and to detect unintended mutations that may occur during strain construction and genome editing.

Recursive genomewide recombination and sequencing reveals a key refinement step in the evolution of a metabolic innovation in Escherichia coli.

Evolutionary innovations often arise from complex genetic and ecological interactions, which can make it challenging to understand retrospectively how a novel trait arose. In a long-term experiment, Escherichia coli gained the ability to use abundant citrate (Cit(+)) in the growth medium after ∼31,500 generations of evolution. Exploiting this previously untapped resource was highly beneficial: later Cit(+) variants achieve a much higher population density in this environment. All Cit(+) individuals share a mutation that activates aerobic expression of the citT citrate transporter, but this mutation confers only an extremely weak Cit(+) phenotype on its own. To determine which of the other >70 mutations in early Cit(+) clones were needed to take full advantage of citrate, we developed a recursive genomewide recombination and sequencing method (REGRES) and performed genetic backcrosses to purge mutations not required for Cit(+) from an evolved strain. We discovered a mutation that increased expression of the dctA C4-dicarboxylate transporter greatly enhanced the Cit(+) phenotype after it evolved. Surprisingly, strains containing just the citT and dctA mutations fully use citrate, indicating that earlier mutations thought to have potentiated the initial evolution of Cit(+) are not required for expression of the refined version of this trait. Instead, this metabolic innovation may be contingent on a genetic background, and possibly ecological context, that enabled citT mutants to persist among competitors long enough to obtain dctA or equivalent mutations that conferred an overwhelming advantage. More generally, refinement of an emergent trait from a rudimentary form may be crucial to its evolutionary success.

Detecting rare structural variation in evolving microbial populations from new sequence junctions using breseq.

New mutations leading to structural variation (SV) in genomes-in the form of mobile element insertions, large deletions, gene duplications, and other chromosomal rearrangements-can play a key role in microbial evolution. Yet, SV is considerably more difficult to predict from short-read genome resequencing data than single-nucleotide substitutions and indels (SN), so it is not yet routinely identified in studies that profile population-level genetic diversity over time in evolution experiments. We implemented an algorithm for detecting polymorphic SV as part of the breseq computational pipeline. This procedure examines split-read alignments, in which the two ends of a single sequencing read match disjoint locations in the reference genome, in order to detect structural variants and estimate their frequencies within a sample. We tested our algorithm using simulated Escherichia coli data and then applied it to 500- and 1000-generation population samples from the Lenski E. coli long-term evolution experiment (LTEE). Knowledge of genes that are targets of selection in the LTEE and mutations present in previously analyzed clonal isolates allowed us to evaluate the accuracy of our procedure. Overall, SV accounted for ~25% of the genetic diversity found in these samples. By profiling rare SV, we were able to identify many cases where alternative mutations in key genes transiently competed within a single population. We also found, unexpectedly, that mutations in two genes that rose to prominence at these early time points always went extinct in the long term. Because it is not limited by the base-calling error rate of the sequencing technology, our approach for identifying rare SV in whole-population samples may have a lower detection limit than similar predictions of SNs in these data sets. We anticipate that this functionality of breseq will be useful for providing a more complete picture of genome dynamics during evolution experiments with haploid microorganisms.

ChIP-seq defined genome-wide map of TGFß/SMAD4 targets: implications with clinical outcome of ovarian cancer.

Deregulation of the transforming growth factor-ß (TGFß) signaling pathway in epithelial ovarian cancer has been reported, but the precise mechanism underlying disrupted TGFß signaling in the disease remains unclear. We performed chromatin immunoprecipitation followed by sequencing (ChIP-seq) to investigate genome-wide screening of TGFß-induced SMAD4 binding in epithelial ovarian cancer. Following TGFß stimulation of the A2780 epithelial ovarian cancer cell line, we identified 2,362 SMAD4 binding loci and 318 differentially expressed SMAD4 target genes. Comprehensive examination of SMAD4-bound loci, revealed four distinct binding patterns: 1) Basal; 2) Shift; 3) Stimulated Only; 4) Unstimulated Only. TGFß stimulated SMAD4-bound loci were primarily classified as either Stimulated only (74%) or Shift (25%), indicating that TGFß-stimulation alters SMAD4 binding patterns in epithelial ovarian cancer cells. Furthermore, based on gene regulatory network analysis, we determined that the TGFß-induced, SMAD4-dependent regulatory network was strikingly different in ovarian cancer compared to normal cells. Importantly, the TGFß/SMAD4 target genes identified in the A2780 epithelial ovarian cancer cell line were predictive of patient survival, based on in silico mining of publically available patient data bases. In conclusion, our data highlight the utility of next generation sequencing technology to identify genome-wide SMAD4 target genes in epithelial ovarian cancer and link aberrant TGFß/SMAD signaling to ovarian tumorigenesis. Furthermore, the identified SMAD4 binding loci, combined with gene expression profiling and in silico data mining of patient cohorts, may provide a powerful approach to determine potential gene signatures with biological and future translational research in ovarian and other cancers.

Aberrant TGFß/SMAD4 signaling contributes to epigenetic silencing of a putative tumor suppressor, RunX1T1 in ovarian cancer.

Aberrant TGFß signaling pathway may alter the expression of down-stream targets and promotes ovarian carcinogenesis. However, the mechanism of this impairment is not fully understood. Our previous study has identified RunX1T1 as a putative SMAD4 target in an immortalized ovarian surface epithelial cell line, IOSE. In this study, we report that transcription of RunX1T1 was confirmed to be positively regulated by SMAD4 in IOSE cells and epigenetically silenced in a panel of ovarian cancer cell lines by promoter hypermethylation and histone methylation at H3 lysine 9. SMAD4 depletion increased repressive histone modifications of RunX1T1 promoter without affecting promoter methylation in IOSE cells. Epigenetic treatment can restore RunX1T1 expression by reversing its epigenetic status in MCP3 ovarian cancer cells. When transiently treated with a demethylating agent, the expression of RunX1T1 was partially restored in MCP3 cells, but gradual re-silencing through promoter re-methylation was observed after the treatment. Interestingly, SMAD4 knockdown accelerated this re-silencing process, suggesting that normal TGF-beta signaling is essential for the maintenance of RunX1T1 expression. In vivo analysis confirmed that hypermethylation of RunX1T1 was detected in 35.7% (34/95) of ovarian tumors with high clinical stages (P=0.035) and in 83% (5/6) of primary ovarian cancer-initiating cells. Additionally, concurrent methylation of RunX1T1 and another SMAD4 target, FBXO32 which was previously found to be hypermethylated in ovarian cancer was observed in this same sample cohort (P< 0.05). Restoration of RunX1T1 inhibited cancer cell growth. Taken together, dysregulated TGFß/SMAD4 signaling may lead to epigenetic silencing of a putative tumor suppressor, RunX1T1, during ovarian carcinogenesis.

Epigenetic silencing mediated through activated PI3K/AKT signaling in breast cancer.

Trimethylation of histone 3 lysine 27 (H3K27me3) is a critical epigenetic mark for the maintenance of gene silencing. Additional accumulation of DNA methylation in target loci is thought to cooperatively support this epigenetic silencing during tumorigenesis. However, molecular mechanisms underlying the complex interplay between the two marks remain to be explored. Here we show that activation of PI3K/AKT signaling can be a trigger of this epigenetic processing at many downstream target genes. We also find that DNA methylation can be acquired at the same loci in cancer cells, thereby reinforcing permanent repression in those losing the H3K27me3 mark. Because of a link between PI3K/AKT signaling and epigenetic alterations, we conducted epigenetic therapies in conjunction with the signaling-targeted treatment. These combined treatments synergistically relieve gene silencing and suppress cancer cell growth in vitro and in xenografts. The new finding has important implications for improving targeted cancer therapies in the future.

Epigenetic influences of low-dose bisphenol A in primary human breast epithelial cells.

Substantial evidence indicates that exposure to bisphenol A (BPA) during early development may increase breast cancer risk later in life. The changes may persist into puberty and adulthood, suggesting an epigenetic process being imposed in differentiated breast epithelial cells. The molecular mechanisms by which early memory of BPA exposure is imprinted in breast progenitor cells and then passed onto their epithelial progeny are not well understood. The aim of this study was to examine epigenetic changes in breast epithelial cells treated with low-dose BPA. We also investigated the effect of BPA on the ERα signaling pathway and global gene expression profiles. Compared to control cells, nuclear internalization of ERα was observed in epithelial cells preexposed to BPA. We identified 170 genes with similar expression changes in response to BPA. Functional analysis confirms that gene suppression was mediated in part through an ERα-dependent pathway. As a result of exposure to BPA or other estrogen-like chemicals, the expression of lysosomal-associated membrane protein 3 (LAMP3) became epigenetically silenced in breast epithelial cells. Furthermore, increased DNA methylation in the LAMP3 CpG island was this repressive mark preferentially occurred in ERα-positive breast tumors. These results suggest that the in vitro system developed in our laboratory is a valuable tool for exposure studies of BPA and other xenoestrogens in human cells. Individual and geographical differences may contribute to altered patterns of gene expression and DNA methylation in susceptible loci. Combination of our exposure model with epigenetic analysis and other biochemical assays can give insight into the heritable effect of low-dose BPA in human cells.

Estrogen-mediated epigenetic repression of large chromosomal regions through DNA looping.

The current concept of epigenetic repression is based on one repressor unit corresponding to one silent gene. This notion, however, cannot adequately explain concurrent silencing of multiple loci observed in large chromosome regions. The long-range epigenetic silencing (LRES) can be a frequent occurrence throughout the human genome. To comprehensively characterize the influence of estrogen signaling on LRES, we analyzed transcriptome, methylome, and estrogen receptor alpha (ESR1)-binding datasets from normal breast epithelia and breast cancer cells. This "omics" approach uncovered 11 large repressive zones (range, 0.35 approximately 5.98 megabases), including a 14-gene cluster located on 16p11.2. In normal cells, estrogen signaling induced transient formation of multiple DNA loops in the 16p11.2 region by bringing 14 distant loci to focal ESR1-docking sites for coordinate repression. However, the plasticity of this free DNA movement was reduced in breast cancer cells. Together with the acquisition of DNA methylation and repressive chromatin modifications at the 16p11.2 loci, an inflexible DNA scaffold may be a novel determinant used by breast cancer cells to reinforce estrogen-mediated repression.

Promoter hypermethylation of FBXO32, a novel TGF-beta/SMAD4 target gene and tumor suppressor, is associated with poor prognosis in human ovarian cancer.

Resistance to TGF-beta is frequently observed in ovarian cancer, and disrupted TGF-beta/SMAD4 signaling results in the aberrant expression of downstream target genes in the disease. Our previous study showed that ADAM19, a SMAD4 target gene, is downregulated through epigenetic mechanisms in ovarian cancer with aberrant TGF-beta/SMAD4 signaling. In this study, we investigated the mechanism of downregulation of FBXO32, another SMAD4 target gene, and the clinical significance of the loss of FBXO32 expression in ovarian cancer. Expression of FBXO32 was observed in the normal ovarian surface epithelium, but not in ovarian cancer cell lines. FBXO32 methylation was observed in ovarian cancer cell lines displaying constitutive TGF-beta/SMAD4 signaling, and epigenetic drug treatment restored FBXO32 expression in ovarian cancer cell lines regardless of FBXO32 methylation status, suggesting that epigenetic regulation of this gene in ovarian cancer may be a common event. In advanced-stage ovarian tumors, a significant (29.3%; P<0.05) methylation frequency of FBXO32 was observed and the association between FBXO32 methylation and shorter progression-free survival was significant, as determined by both Kaplan-Meier analysis (P<0.05) and multivariate Cox regression analysis (hazard ratio: 1.003, P<0.05). Reexpression of FBXO32 markedly reduced proliferation of a platinum-resistant ovarian cancer cell line both in vitro and in vivo, due to increased apoptosis of the cells, and resensitized ovarian cancer cells to cisplatin. In conclusion, the novel tumor suppressor FBXO32 is epigenetically silenced in ovarian cancer cell lines with disrupted TGF-beta/SMAD4 signaling, and FBXO32 methylation status predicts survival in patients with ovarian cancer.

Epigenetic repression of microRNA-129-2 leads to overexpression of SOX4 oncogene in endometrial cancer.

Genetic amplification, mutation, and translocation are known to play a causal role in the upregulation of an oncogene in cancer cells. Here, we report an emerging role of microRNA, the epigenetic deregulation of which may also lead to this oncogenic activation. SOX4, an oncogene belonging to the SRY-related high mobility group box family, was found to be overexpressed (P < 0.005) in endometrial tumors (n = 74) compared with uninvolved controls (n = 20). This gene is computationally predicted to be the target of a microRNA, miR-129-2. When compared with the matched endometria, the expression of miR-129-2 was lost in 27 of 31 primary endometrial tumors that also showed a concomitant gain of SOX4 expression (P < 0.001). This inverse relationship is associated with hypermethylation of the miR-129-2 CpG island, which was observed in endometrial cancer cell lines (n = 6) and 68% of 117 endometrioid endometrial tumors analyzed. Reactivation of miR-129-2 in cancer cells by pharmacologic induction of histone acetylation and DNA demethylation resulted in decreased SOX4 expression. In addition, restoration of miR-129-2 by cell transfection led to decreased SOX4 expression and reduced proliferation of cancer cells. Further analysis found a significant correlation of hypermethylated miR-129-2 with microsatellite instability and MLH1 methylation status (P < 0.001) and poor overall survival (P < 0.039) in patients. Therefore, these results imply that the aberrant expression of SOX4 is, in part, caused by epigenetic repression of miR-129-2 in endometrial cancer. Unlike the notion that promoter hypomethylation may upregulate an oncogene, we present a new paradigm in which hypermethylation-mediated silencing of a microRNA derepresses its oncogenic target in cancer cells.

Methylation analysis by microarray.

Differential methylation hybridization (DMH) is a high-throughput DNA methylation screening tool that utilizes methylation-sensitive restriction enzymes to profile methylated fragments by hybridizing them to a CpG island microarray. This array contains probes spanning all the 27,800 islands annotated in the UCSC Genome Browser. Herein we describe a DMH protocol with clearly identified quality control points. In this manner, samples that are unlikely to provide good read-outs for differential methylation profiles between the test and the control samples will be identified and repeated with appropriate modifications. The step-by-step laboratory DMH protocol is described. In addition, we provide descriptions regarding DMH data analysis, including image quantification, background correction, and statistical procedures for both exploratory analysis and more formal inferences. Issues regarding quality control are addressed as well.

Xenoestrogen-induced epigenetic repression of microRNA-9-3 in breast epithelial cells.

Early exposure to xenoestrogens may predispose to breast cancer risk later in adult life. It is likely that long-lived, self-regenerating epithelial progenitor cells are more susceptible to these exposure injuries over time and transmit the injured memory through epigenetic mechanisms to their differentiated progeny. Here, we used progenitor-containing mammospheres as an in vitro exposure model to study this epigenetic effect. Expression profiling identified that, relative to control cells, 9.1% of microRNAs (82 of 898 loci) were altered in epithelial progeny derived from mammospheres exposed to a synthetic estrogen, diethylstilbestrol. Repressive chromatin marks, trimethyl Lys27 of histone H3 (H3K27me3) and dimethyl Lys9 of histone H3 (H3K9me2), were found at a down-regulated locus, miR-9-3, in epithelial cells preexposed to diethylstilbestrol. This was accompanied by recruitment of DNA methyltransferase 1 that caused an aberrant increase in DNA methylation of its promoter CpG island in mammosphere-derived epithelial cells on diethylstilbestrol preexposure. Functional analyses suggest that miR-9-3 plays a role in the p53-related apoptotic pathway. Epigenetic silencing of this gene, therefore, reduces this cellular function and promotes the proliferation of breast cancer cells. Promoter hypermethylation of this microRNA may be a hallmark for early breast cancer development, and restoration of its expression by epigenetic and microRNA-based therapies is another viable option for future treatment of this disease.